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In quantum spin chains at criticality, two types of scaling for the entanglement entropy exist: one 
comes from conformal field theory (CFT), and the other is for entanglement support of matrix prod- 
uct state (MPS) approximation. They indicates that the matrix dimension of the MPS represents a 
length scale of spin correlation. On the other hand, the quantum spin-chain models can be mapped 
onto two-dimensional (2D) classical ones. Motivated by the scaling and the mapping, we introduce 
new entanglement entropy for 2D classical spin configuration as well as entanglement spectrum, and 
examine their basic properties in Ising and 3-state Potts models on the square lattice. They are 
defined by the singular values of the reduced density matrix for a Monte Carlo snapshot. We find 
scaling relations concerned with length scales in the snapshot at Tc. There, the spin configuration 
is fractal, and various sizes of ordered clusters coexist. Then, the singular values automatically 
decompose the original snapshot into a set of images with different length scale. This is the origin 
of the scaling. In contrast to the MPS scaling, long-range spin correlation can be described by 
only few singular values. Furthermore, we find multiple gaps in the entanglement spectrum, and in 
contrast to standard topological phases, the low- lying entanglement levels below the gap represent 
spontaneous symmetry breaking. Based on these observations, we discuss about the amount of 
information contained in one snapshot in a viewpoint of the CFT scaling. 

PACS numbers: 05.10.Cc, 89.70.Cf, 11.25.Hf 



I. INTRODUCTION 

The entanglement entropy is a common language con- 
necting among various fields such as quantum informa- 
tion, quantum gravity, and condensed matter physics. 
The main reason of this wide applicability comes from 
a fact that the entropy picks up universality irrespec- 
tive of details of their models. The entropy represents 
the amount of information across the boundary between 
a subsystem A of linear size L and its environment B. 
Starting with a wave function of the total system |?/^), we 
first define the density matrix of A by pA = tr^ IV^) (V^l, 
which traces out degree of freedom inside of B. Then, 
the entanglement entropy Sa is given by 



Sa = -tr^ (pa log Pa)- 



(1) 



It has been extensively examined how this entropy be- 
haves as functions of L and spatial dimension d. 

A well-known formula is called 'area-law scaling', S ex 
which tells us non-extensivity of S in contrast to 
the thermal entropy. This formula was originally in- 
troduced in a context of black-hole physics (Bekenstein- 
Hawking entropy) [i|-[3|, and examination of the scaling 
and its violation has been a hot topic in condensed mat- 
ter physics [4-11]. 

The violation occurs in cases of one-dimensional (ID) 
critical systems and models with Fermi surface. In these 
cases, the scaling contains logarithmic correction. 
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(2) 



where C is related to the number of excitation modes 
across the boundary between A and 5, and is equal to 



the central charge c of conformal field theory (CFT) in 
d = 1. In the CFT, the entropy is roughly given by a 
logarithm of a two-point correlation function for scaling 
operators, and thus Eq. ([2j) naturally appears. Away 
from a critical point, the entropy is deformed as 



S=-cA\og(^, 
6 



(3) 



with correlation length ^ and the number of boundary 
points A of A. 

On the other hand, there is another type of entropy 
scaling which does not contain the universality parame- 
ter c and any length scale explicitely. When we take ID 
quantum critical models by using variational optimiza- 
tion of matrix product state (MPS) [12], it is conjectured 
that the half-chain entanglement entropy behaves as 
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logx, 



(4) 



where x is matrix dimension of MPS. This conjecture 
was recently found in two specific models with different 
central charges, respectively: one is transverse-field Ising 
chain (c = 1/2), and the other is 5* = 1 XX Z chain with 
uniaxial anisotropy (c = 1) [isl, ll^. The logx depen- 
dence on this entropy can be interpleted as a result of 
quantum entanglement between A and B. This is be- 
cause = logx maximally entangled-pair state, 
lip) — (1/y^) X]n=i l^^)^ where the states in A and B 
are labeled by n and n, respectively, and one of the two 
degrees of freedom, n or n, is traced out. Here, \\\)) is 
also a particular form of MPS. Furthermore, a prefactor 
1 /6 is expected to be a character of the Virasoro algebra 
in CFT, although the microscopic understanding has not 
been obtained yet. 
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In general, MPS for x = 1 represents local approxi- 
mation (no entanglement, =0), and is asymptotically 
exact if we could take a sufficiently large x value. Then, 
approximately taking a finite x value would limit spin 
correlation or quantum entanglement to finite-range one. 
In that sense, x controls the range of the spin correla- 
tion. By combining Eq. (jH with Eq. ([3j), we know that 
the effective correlation length of MPS is given by 

Ceff = x'/^ (5) 

where we take A=l because Eq. (|4]) is obtained for the 
half of an infinitely-long chain. Therefore, this x value is 
actually related to the length scale that represents how 
presicely the spin correlation is taken into account. How- 
ever, this is somehow strange, since x is just a parameter 
for how many singular values of the matrices in MPS are 
taken. The above consideration suggests that the length- 
scale control given by Eq. (|5]) is a fundamental function 
of the singular value decomposition (SVD). The issue to 
be resolved here is why the SVD automatically produces 
the length scale. 

Here, we address this issue in a viewpoint of quantum 
/ classical correspondence. Usually, a ID quantum spin 
model is transformed into a 2D classical spin model by 
the Suzuki- Trotter decomposition. Then, we can handle 
a Monte Carlo (MC) simulation, and obtain a snapshot of 
particular spin configuration. The correspondence may 
predict that a length scale characterized by x the 
quantum side is hidden in the snapshot. Therefore, we 
attempt to search the hidden length scale, and discuss 
about physical meaning of Eq. (|5|). This is a purpose of 
this paper. 

For this purpose, we introduce new entanglement en- 
tropy for a snapshot calculated by a MC simulation. 
This is the von Neumann entropy defined by the singular 
values of the reduced density matrix for the snapshot. 
Then, we find two scaling relations of the entropy in the 
Ising and the 3-state Potts models that are analogous to 
Eq. Q (or Eq. (|2|)) and Eq. (g]). Furthermore, the scal- 
ing also appears on the positions of the multiple gaps in 
the entanglement spectrum. A key factor for the scaling 
is fractal spin configuration at Tc. The two scaling rela- 
tions come from short- and long-range spin correlation in 
the fractal. A role of the SVD on the length-scale con- 
trol is to decompose the original snapshot into a set of 
images with different length scales, respectively. Then, 
each scale is characterized by one of the multiple gaps in 
the entanglement spectrum. We discuss about similar- 
ity and difference between the new entropy scaling and 
standard one in ID quantum systems, and also discuss 
about possible presence of a topological term hidden in 
our scaling relation. 

The paper is organized as follows. In Sec. II, we de- 
fine the entanglement entropy for a snapshot, which is 
a key ingredient in this paper, and present outline of 
our method. Then, in Sec. HI, basic properties of the 
entanglment entropy and the entanglement spectrum for 
square-lattice Ising ferromagnet are presented. The main 



objective is to show temperature and system-size depen- 
dence of the entropy as well as the spectrum in order 
to extract scaling relations. We also examine entangle- 
ment support of our method by changing the number 
of the singular values which are taken into account. In 
Sec. IV, coarse-grained snapshots are shown, and we dis- 
cuss about the key mechanism of the length-scale control 
hidden in the SVD. In Sec. V, we examine the 3-state 
Potts model in order to conffim universality of our scal- 
ing relations obtained in the analysis of the Ising model. 
In Sec. VI, we discuss about the topological entanglement 
entropy in a viewpoint of the entanglement gap. Finally, 
we summarize our results. 



II. METHOD 

We start with the Ising model on the square lattice: 
H = -J^a,aj. (6) 

where cr^ = ±1 and the sum runs over the nearest neigh- 
bor lattice sites (i, j), and J(> 0) is exchange interaction. 
The system size is taken to be L x L. According to dual 
transformation, the critical temperature is known to be 
Tc/J = 2/log(l + V2) = 2.2692. The central charge of 
the Ising model is c = 1/2. 

First, we are going to obtain a snapshot of a spin con- 
figuration m{x,y) = <Ji with i = {x^y). We can freely 
choose a method for obtaining the snapshot. Here, we 
will use MC simulation. We regard m(x, y) as a matrix, 
and calculate the reduced density matrices defined by 

px{x,x') = ^m{x,y)m{x\y), (7) 

y 

pY{y,y) = ^m{x,y)m{x,y), (8) 

X 

where we trace over y (x)-component in px {px)- Let us 
decompose m{x^y) into a set of the sigular values {A^} 
and the column unitary matrices {Un{x)} and {Vn{y)}'. 

L 

m{x,y) = '^Un{x)y/A^Vn{y). (9) 

n=l 

Mathematically, {A^} are uniquely determined, while 
{Un} and {Vn} are not. Thus if some universal features 
could be extracted from a snapshot, those should be rep- 
resented by a function of A^. By substituting Eq. ([9]) 
into Eqs. ([7]) and (j8j), we have 

L 

px{x,x') = Y,Un{x)KnUn{x), (10) 
L 

pY{y,y') = Y.^n{y)KVn{y'). (ii) 

n=l 
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Thus, the set of {A^} is obtained by diagonahzing px or 
pY- It is noted that the eigenvalue spectrum of px is the 
same as that of py. Even if we consider a rectangular 
lattice with M x N sites, the nonzero eigenvalues of px 
and pY are the same, and the number the eigenvalues 
is L = min(M, A^). We align the eigenvalues so that 
Ai>A2>--->Al. Each eigenvalue A^ is normalized 
to be An = An /C with a constant C so as to satisfy 

L 

J2Xn = l. (12) 

n=l 

Next, we define the von Neumann entropy of a snap- 
shot analogous to Eq. ([T]). That is the amount of entan- 
glement between x- and ^-components defined by 

X 

^X = -IZ^nlogAn, (13) 

n=l 

with L. We abbreviate Sl to S. This is a key quan- 
tity in this study. However, at the present stage, we do 
not know whether this is related to the standard entan- 
glement entropy shown in Eqs. (|2]) and Hereafter, we 
will present basic properties of S in detail. Before going 
to the detail, it is theoretically clear that this entropy 
becomes maximum when we take = 1/x for any n. 
Then we have 

S^<-Y,-^^Z-=^ogx. (14) 

and 5*^ is bounded by log x- 

I have performed MC simulation by a standard 
Metropolis algorithm in order to obtain m(x, y). Periodic 
boundary condition is taken into account for the square 
lattice. Starting with the temperature T = 3.02J, 10^ 
MC steps are taken for thermal equillibrium. Here, one 
MC step counts L x L updates. After that, I gradually 
reduce temperature by AT = 0.05J and take 10^ MC 
steps for convergence at each T. When calculating tem- 
perature dependence of the total entropy S in detail, I 
take AT = 0.01 J, and in this case the MC step in each T 
is taken to be 10^ ~ 10^ depending on the system size. I 
have also confirmed numerical convergency by taking 10^ 
MC steps for some T values. I have observed snapshots 
and their entropy across Tc. 

III. ENTANGLEMENT ENTROPY AND 
ENTANGLEMENT SPECTRUM OF THE ISING 
MODEL 

A. Scaling relation for the entropy at T ^ Tc 

Let us look at Fig. [1] where basic properties of the 
entropy S are summarized. Since we would like to ex- 
amine the amount of information in one snapshot, we do 
not take statistical average of S except for cases that we 
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FIG. 1: (a-c) Temperature and system-size dependence on 
the entropy S and a = S — logL: L = 64 = 2^ (blue circles), 
L = 128 = 2^ (purple circles), L = 256 = 2^ (red circles), 
L = 512 = 2^ (black circles), L = 1024 = 2^° (solid line). 
In pannel (c), we have avaraged 10^ samples for L = 64. A 
dashed vertical line is a guide to Tc. (d) finite-size scaling 
for a near Tci T = 2. 25 J (filled triangles), T = 2. 26 J (filled 
circles), T = 2. 27 J - Tc (open circles), T = 2. 28 J (filled 
diamonds), and T = 2. 29 J (crosses). We have avaraged 10^ 
and 10^ samples for L < 256 and L = 512, respectively. Note 
that a relative statistical error between 10^ and 10^ samples 
is Aa - 0.04 for L = 256 at T = 2.29 J. 



need precise scaling. Fortunately, the statistical error of 
S becomes smaller as L increases as shown in Figs. [1] (a) 
and (b), and the effect of self-avarage on S seem to 
be much better than that of thermodynamic quantities. 
This small variance guarantees that our data do not suf- 
fer from the severe statistical error at least for large-T 
region. 

As for Tc, temperature dependence of 5 is a good mea- 
sure, since S behaves quite differently below and above 
Tc as shown in Fig. [T](a). Below Tc, the system is in the 
ferromagnetically ordered state (Ai = 1 and otherwise 
0), and then S should go to zero. Above Tc, the spin 
configuration is paramagnetic, leading to high entropy. 
In that sense, the T-dependence is similar to that of the 
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thermal entropy. We see that S' for L = 1024 largely 
drops at 2.26J < T < 2.27J with decreasing tempera- 
ture, suggesting phase transition. This position is very 
close to the exact Tc. Furthermore, at T > Tc, we find 



S = \oq:L 



(15) 



with a < 0. Figure [T] (b) plots a = S — logL instead 
of S in order to show this scaling clearly. In the next 
paragraph, we will obtain a = — 7r/4 ^ —0.77 in the 
large-T and the large-L Hmits. Actually, a at T = 3.02 J 
is close to — 7r/4. At Tc, the value is estimated to be 
a ~ — 2 as shown in Fig.[T](c), where the data for L = 64 
are avaraged by 10^ samples. We pick up each sample per 
1 MC step after thermalization. The finite-size scaling for 
a near Tc is also presented in Fig. [1] (d), and the result 
supports a ^ —2. 

However, the logL dependence on S at and above Tc 
seems to come from different origins. In Fig.[T](d), the a 
value above Tc (T = 2.29 J) increases slightly with L, and 
finally converges for L ^ 512. On the other hand, the a 
value gradually decreases with increasing L at T = 2.27 J. 
Furthermore, the eigenvalue distribution has strong T- 
dependence particularly near Tc as shown later. In the 
following subsections, we examine the physical origins of 
the scaling Eq. ([15]) at and above Tc separately. 



B. Random matrix theory in the large-T limit 



S is precisely determined by random matrix theory, since 
the spin configuration is paramagnetic (random). Here, 
we introduce the eigenvalue distribution 

1 ^ 

n=l 

1^1 

7^0+ T ^ TT (X - An)2 +7^ 

n=l ^ ^ 

and according to the random matrix theory p{x) should 
asymptotically approach the Mar cenko- Past ur law in the 
large-T limit 



1 



p(,) = ____y^(4T^, 



(18) 



for < X < 4T with variance T. Actually, the numeri- 
cally obtained distribution for T = 10 J and T = 512 well 
fit with this equation as shown in Fig. [2] (b). For those 
parameters, the average of off-diagonal components of 
px is 0.2662 which is very small, and the variance of the 
off-diagonal components is 530.07 ~ T. The maximum 
eigenvalue is A^«^ = 2078.125 - 4T. These data also fit 
with the random matrix theory. Then, S can be evalu- 
ated as follows: We transform Eq. ([13]) with x = ^ i^i^o 
an integral form with use of Eq. (p!6|) 
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FIG. 2: (a) Normalized eigenvalues An (L = 512): T = 10.0 J 
(black) and T = 40 J (red circles), (b) Eigenvalue distribution 
p{x) {L = 512): T = 10 J (black circles). A red solid line rep- 
resents the asymptotic distribution curve by the Marcenko- 
Pastur law. A broadening factor in Eq. (pT|) is taken to be 
= 10 < /^rnax ^ rpj^^ -^^^^ shows the numerically ob- 
tained entropy 5* as a function of log L at T = 10 J. The MC 
steps are lO'^. A guide line in the inset represents Eq. (|20)) . 

Let us first examine the eigenvalues above Tc. In 
Fig. [2] (a), the eigenvalue as a function of n decays 
slowly, and then all of the eigenvalues play a role on the 
entropy. This behavior is unchanged for large-T region, 
and we see that the data with T = 10 J and 40 J are al- 
most the same. In the large-T limit, the upper bound of 



^ dx^^(5(x-An)glog(^) 
dxLp{x)^ log (^^^ , 



(19) 



and substituting Eq. (p!8|) to this equation. Then, we ob- 
tain the following scaling relation: 



S 



^0 



V^(4T - x) 



dx- 



27rT2 



log 



loe;T 



log (2(1 + t)) 



TT 



where C is taken to satisfy 

L f.AL 



C = A„ = / dxp{x)x = L^. 



(20) 



(21) 



We have confirmed that this relation is strictly hold in our 
numerical simulation. Therefore, the log L dependence in 
the large-T limit comes from the universal feature of the 
random matrix, and the diviation from the completely 
random state is characterized by a. Numerical data in 
the inset of Fig. [2] (b) support this analytical result. 
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C. Eigenvalue distribution and Entanglement 
spectrum at Tc 




n n 



FIG. 3: (a) Normalized eigenvalues An (L = 512): T = 2. 27 J 
(blue), T = 3.02 J (red), and T = 10.0 J (black). The in- 
set is An for small n region, (b) Eigenvalue distribution 
p{x) (L = 512): T = 2. 27 J (black open circles). A red 
solid line represents the asymptotic distribution curve by the 
Marcenko-Pastur law. We take 7 = 10. (c) Entanglement 
spectrum for T = 2. 27 J: L = 128 (red), L = 256 (purple), 
and L — 512 (black). We plot 10 different-sample data for 
each parameter on the same pannel. (d-f) Logarithmic plot 
of eigenvalue An for L = 128, 256, and 512. 

On the other hand, the nature of the eigenvalues at 
Tc is quite different from that above Tc. We present the 
eigenvalues and the distribution function in Figs. [3] (a) 
and (b). In Fig. [3] (a), we observe evolution of A^ in smah 
n region toward Tc, and finally Ai ^ 1 at zero tempera- 
ture. Thus, the data at Tc can be viewed as a mixture of 
various length scales, and small and large n regions rep- 
resent long-range (ferromagnetic) and short-range (para- 
magnetic) components of spin correlation, respectively. 
As shown in Fig. [3](b), the eigenvalue distribution devi- 
ates from that of the random matrix, and thus the log L 
dependence at Tc is not due to the random matrix. Here, 
the average and the variance of off-diagonal components 
of px are 92.022 > and 5612.122 > L, respectively, 
and A^^^ = 53350.01 > 4L. These data also differ from 



those expected by the random matrix theory. 

The separation of A^ into different length scales is more 
clearly seen in the entanglement spectrum [l5| defined by 

En = - log An. (22) 

In Fig. [3] (c) and the inset, we plot 10 different-sample 
data for each parameter in order to reduce finite-size ef- 
fects. We find multiple entanglement gaps near Tc 0. 
In addition to large gaps at around En = 2 and 3, gap-like 
anomalies also appear at En ^ 3.75 and 4.25 for L = 512. 
Because of the presence of the multiple gaps, the states 
are separated into a set of different entanglement levels. 
The state below the largest gap represents a precursor of 
ferromagnetic long-range order, while nearly continuum 
states far above the large gaps represent paramagnetic 
short-range spin configuration. In between, each sector 
separated by the gaps would have intermediate length 
scales of spin correlation. Then, we expect that the mul- 
tiple gaps characterize various length scales, and their 
mixture leads to critical behavior. 

Let us look at more about the gap- like anomalies above 
En = 3 in Figs. [3] (d-f), where log^^g ^ En is plotted. 
We find a kink structure, and the position of the kink, 
Ukinki exactly traces one of the multiple entanglement 
gaps. The kink position increases as L. We observe that 
^kink = 3, 6, and 9 for L = 2^, 2^, and 2^, respectively. 
Thus, this observation tells us 

\rikink = log2 L - 6, (23) 

for L > 2^. The kink position gives us a criterion of 
how many eigenvalues for the given L play an important 
role on the spin correlation associated with the critical 
behavior. As we increase L, the discreteness of the lat- 
tice tends to disappear, and eventually various sizes of 
ferromagnetic clusters appear at Tc. Therefore, Ukink 
increases as L. The kink structure does not appear in 
Fig. [2] (a) for large T, and thus Eq. ([23]) is a character 
which only appears near Tc. We consider that this logL 
dependence on Ukink is closely related to the logL de- 
pendence on S in Eq. (p!5]) . It should be noted that only 
few states are concerned with the long-range spin corre- 
lation of the critical behavior in our expression. This is 
quite contrast to the MPS scaling, and this point will be 
important in the later discussion. 

D. Compasiron with the thermal entropy 

In order to see the pecuriality of S near Tc, it is mean- 
ingful to compare S with the thermal entropy ^^^ermoZ 
The thermal entropy also goes to zero below Tc due to 
ferromagnetic order, while in the large-T limit we have 
a high constant value log 2 coming from the number 
of all possible spin configurations 2^ . Thus, the physi- 
cal meaning of S and S^^^^^^^ is essentially the same in 
the both limits. Figure 3] shows comparison between S 
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thermodynamic limit, s^^^™^^ = Hirl^oo S^^^^^^^ / L'^ ^ 
is given by the Onsager's solution. First, we transform 
the thermodynamic first law into the following form 



Aherraal 



(24) 



with free energy per site / and inverse temperature P = 
1/T. Then we substitute the Onsager's free energy into 
the above equation to obtain ^^^ermaZ. 

~l3f = \ log(2 sinh 2K) + dqe{q, K), (25) 

where K = /3J, and e(g, K) is a solution of the following 
equation 



cosh e(g, K) = cosh 2K coth 2K — cos q. 



(26) 




FIG. 4: Comparison between 5'/(log L — 7r/4) (solid line, L — 
1024) with s*^"""'"7log2 (dashed line, exact). 

In Fig. m we show the normalized data so that both of 
them approach unity in the large-T limit. We find that 
the change in S near Tc is more cusp-like in comparison 
with ^^^erma/ xhis cusp-likc feature is an evidence that 
separates the logL dependence near Tc from that above 
Tc. The cusp reminds us a fact that the entanglement 
entropy in transverse-field Ising model increases toward 
the quantum critical point. 



E. Finite-x scaling 

The scaling Eq. ([15]) at Tc can be understood by calcu- 
lating X dependence on 5*^. Figures [5] and [6] show snap- 
shots of spin configurations and their entropy 5*^ at var- 
ious T. We take L = 1024, since we know that a is 
almost converged for this size as discussed in Fig. [H We 
find that the entropy asymptotically approaches 



(27) 



for X ^ '^kink as we increase L. The residual entropy 7' 
does not depend on L, and 7' = 6'i ~ 0.35. Now, we have 
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FIG. 5: Snapshots of spin configurations (L — 256): (a) 
T/J = 3.02, (b) T/J = 2.32, (c) T/J = 2.27, (d) T/J = 1.52. 
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FIG. 6: (a) for L = 1024 at T = 2.27J - Tc (purple 
circles). For compasiron, we also plot the data for L = 128 
(purple box) and L = 512 (purple triangles) at T = 2.21 J . 
A solid line with slope 1/6 is a guide to the eye. A dashed 
line represents — logx — 2. The inset shows ^S'^/^logx 
for L — 1024. The black filled circle represents 1/6. (b) 
Temperature dependence of {L = 1024): T = 3.02 J (red) 
and T = 1.52 J (blue). 



nkink = 12 and logUkink = 2.485 for L = 1024. In the 
inset of Fig. [6] (a), we plot the first derivative of 5*^ with 
respect to logx order to see the slope more precisely. 
Furthermore, for x > ^kinki we find 



Sx = logx + 7, 



(28) 
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with 7 = -2 for L = 1024. We rewrite Eq. (j28j) as 

Sx = l log nfeinfe + 7 + A5^. (29) 

As already mentioned, the additional term A^^ comes 
from short-range spin correlation. Since this component 
originates from paramagnetic spin configuration above 
Tc, the absolute value is larger than the (1/6) logx term. 

Let us first look at Eq. ([271) in Fig. El (b). In the 
ferromagnetically-ordered phase, the entropy should dis- 
appear, and actually the slope of dit T = 1.52 J < Tc 
is almost zero. As we approach T = 2.27 J ^ Tc^ the 
slope of 5*^ for x ^ ^kink gradually increases toward 1/6. 
However, the slope starts to decrease beyond Tc. Actu- 
ally, the slope becomes very small for T = 3. 02 J > Tc. 
According to the previous subsection, the slope near Tc 
would be due to the long-range spin correlation. 

Next we consider Eq. ([28l) . We have already observed 

7 = -2 = a, (30) 

and then the scaling (|28l) after taking x = L is considered 
to be the origin of Eq. (fT5l) at Tc. In Fig[6l somehow 
saturates at around x = but this is finite-size effect. 
We have numerically confirmed that the saturation re- 
duces as we increase L. Thus, for x Eq. (|28l) is 
strictly satisfied in the large-L limit. This feature is also 
consistent with the previous result that S asymptotically 
approaches Eq. ^ib^ in the large-L limit at Tc. With in- 
creasing T, is still proportional to logx foi" X ^ ^kink^ 
while the slope gradually increases. At T = 3.02 J > Tc, 
the slope is about two times larger. Since the residual 
term 7 is a large negative value, finally approaches 
log!/ -\- a dit X = L. 

All of the results presented in this section suggest that 
the origins of the log L dependence in Eq. are clearly 
different at and above Tc. 

F. Block-spin transformation and scaling 

In order to examine the physical meaning of Eqs. ([271 
[3Ql) . it is efficient to introduce the block-spin transforma- 
tion and calculate the entanglement entropy associ- 
ated with the transformation. The block-spin transfor- 
mation merges 3x3 lattice sites together, and the effec- 
tive Ising spin on the new site (the number of the new 
sites is represented by L*) also takes 1 or —1 depending 
on spin configuration that more than half are up or down 
spins, respectively. We continue the transformation un- 
til the system becomes one effective site. Then, /S* for 
L* ^ 1 represents the entropy of the fixed point of this 
renormalization group. The entropy 5* is calculated by 

L* 

5* = -^ a: log a:. (31) 

n=l 

where A* is the eigenvalue of the reduced density matrix 
in the coarse-grained system with the system size L*. 




logL 



FIG. 7: Entanglement entropy 5** for snapshots generated by 
the block-spin transformation at T = 3. 02 J (red), T = 2. 27 J 
(purple), and T = 1.52 J (blue). The original snapshot is 
taken for L — 243 = 3^. We avarage lO'^ samples for the 
calculation of S**. Solid and dashed guide lines are 5** = 
logL* - 2 and 6'* = (1/6) logL*, respectively 

Figure [71 shows L*-dependence on 6** at various T . At 
Tc, we again find the scaling 

S* r.- \o^L\ (32) 
6 

for small L* region, and 

5* -logL* -2, (33) 

for large-L* region. A crossing point of these lines is not 
far from Ukink- Thus, we see that L* plays a similar 
role on x in Eqs. ([27[) and ([28]) . Therefore, it is rea- 
sonable to consider that x represents a length scale as- 
sociated with the couarse graining. We also understand 
that the spin fluctuation in short-range scale has been al- 
ready renormalized by the block-spin transformation in 
small L* region, and only large-scale phenomena survive. 
Then, Eq. ([27|) characterizes this large-scale phenomena. 

IV. COARSE-GRAINED SNAPSHOT AND 
LENGTH SCALE 

A. Coarse-grained snapshot 

In order to understand the nature of the scaling re- 
lations in detail, it is important to remenber that fer- 
romagnetic islands in the snapshot are fractal-like at Tc 
(see Fig. [51(c)). Actually, the system is self-similar, and 
we can always observe various sizes of the islands even 
when we continue to zoom in the system. The variety of 
the island sizes plays a crucial role on the presence of the 
scaling relations. This should be also related to the pres- 
ence of the multiple gaps in the entanglement spectrum. 
However, let us also remember the original definition of 
the entropy 5*^, where x is just the truncation number 
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of the SVD and does not seem to connect directly to any 
length scale. Thus, we should examine how such a length 
scale is automatically generated by the SVD. 

We again call SVD of the snapshot m{x,y) = 

Ylt=i Un{x)\^X^Vn{y). Our target is to look at a 'de- 
formed' snapshot which is defined by restricting the sum 
upto X ill the SVD: 

X 

m^{x,y) = Un {x) \fKiVn {y) , (34) 

n=l 

where m{x^y) = mL{x^y). We expect that this should 
be a 'coarse-grained' image of the original snapshot, if x 
characterizes a length scale. 

Figure [8] (a) is a target snapshot at T = 2.27 J ~ Tc, 
and Figs. [8] (b)-(h) zoom in a region 25 < x < 75 and 
206 <y< 256 of m^{x^y) for various x- First, compar- 
ing Fig. [8] (a) with (b), we see that the original snapshot 
is fractal-like: even if we zoom in the snapshot, various 
sizes of ferromagnetic islands appear sequentially. De- 
creasing the value of we find that global structures 
do not change so much, but y) gradually loses fine 

structures. The snapshot contains various sizes of the 
islands near Tc, and thus the lost of the fine structures 
means that larger and larger islands are damaged by the 
reduction of x- Let us look at one of the smallest islands 
(surrounded by a pentagon) that disappears in pannel 
(e). We also concentrate on the island surrounded by a 
circle. This island is a little bit larger than the smallest 
one. The island still remains in pannel (e), although the 
shape is deformed. Therefore, x is actually controlling 
the accesible length scale of our model. 

B. Layered structure with different length scales 

As already discussed, the length scale is characterized 
by the kink structure (multiple gap structure) in the en- 
tanglement spectrum. In order to see the effect of the 
structure on the snapshot, we present contour map for 
each layer of m(x, y) defined by 

m^''\x,y) = U„{x)^/^V„{y), (35) 

where we call the label n as 'layer'. The parameters are 
taken to be L = 256 and T = 2.27J. In Fig. [9] (b) with 
n = 2, we find that m^'^\x^y) ~ 1 in wide {x^y) region 
(see the vertical axis). When we compare Fig.[9](b) with 
(a), it is clear that the region represents the largest fer- 
romagnetic islands. In Fig. [9] (c) with n = 4, we observe 
many broad peaks, and the positions of some of the repre- 
sentative peaks are assigned to be those of relatively large 
ferromagnetic islands. In Fig.[9](d) with n = 8, the peaks 
become very sharp, and they represents small ferromag- 
netic islands. Further increasing the number of n, we can 
observe much shaper structures representing smaller and 
smaller islands. Since the spatial distribution of a single 
peak corresponds to the size of a ferromagnetic island. 




FIG. 8: Zooming in a region 25 < x < 75 and 206 < y < 256 
of y) for T = 2. 27 J and L = 256: (a) original snapshot 

(L = 256), (b) X = L, This image corresponds to white- 
square region of (a), (c) x = ^/2, (d) x = (e) x = 
(f) x = (g) X = L/32, and (h) x = ^/64. 



we understand that each layer has its own length scale. 
This is the reason why x is a scaling parameter. 

Let us remember that the kink structure in the entan- 
glement spectrum appears at rikink = 6 for L = 256 = 2^. 
Then, we actually see that the data with n = 8 > Ukink 
only represent local components. In the present case, a 
small n represents large ferromagnetic islands, while in 
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FIG. 9: Contour map of m}^\x^y)\ L — 256, rikink — 6, 
and T = 2.27 J ^ (a) original snapshot, (b) n = 2 < rikink^ (c) 
n = 4 < Ukink, and (d) n = 8 > n/einfc- 



FIG. 10: Contour map of Trv^\x^ y)\ L — 256, Ukink — 6, and 
T = 2.27 J ^ (a) original snapshot which differs from Fig. [9] (a), 
(b) n = 2 < Ukink, (c) n = 4 < rifcinfc, and (d) n = 8 > Ukink- 
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Eq. (j5j) a large x value represents successive treatment 
of long-range correlation. Thus, the role of x on ^ is re- 
versed between quantum and classical sides. Figure [TOl 
shows contour map for a snapshot with rather different 
island structures. The tendency of length-scale control 
also appear in this case. When we compare Fig. [9] with 
Fig. [TOl the average size of each island for a fixed n value 
is similar in both Figures. This clearly shows the pres- 
ence of the multiple entanglement gaps. 



Intuitive understanding of the length-scale 
control by SVD 



Large island 



Small island 




Sign 
cliange 



FIG. 11: 2x4 lattice model for coexistence of small and 
large islands: cr{x,y) = +1 (red pixel) and (j{x,y) = —1 
(blue pixel). The left half of the system represents a large 
ferromagnetic island, while m(l,4) represents a small island. 
They are separated by background spins with a{x,y) = —1. 
In this case, the area of the background spins are larger than 
that of the small island. 

Let us explain the key mechanism of the length scale 
control hidden in SVD. Here, we introduce a snapshot 
with 2x4 pixels shown in Fig. [TTJ This is a minimal 
model of the snapshot where small and large ferromag- 
netic islands coexist. The magnetic moment m of this 
snapshot is represented by the matrix 



11-11 
11-1-1 



(36) 



Here, we call m(l,4) as small ferromagnetic island, and 
call the left half as large ferromagnetic island. Then, the 
reduced density matrices px and pY are given by 



Px = 



PY 




(37) 
(38) 



Their nonzero eigenvalues are the same, and they are 
Ai = 6 and A2 = 2. In the present ideal case, the density 
matrix py can be completely decoupled into two sub- 
spaces: /)y(4,4) represents the small island and m(2,4), 
while the remaining 3x3 submatrix represents the large 



island and two background spins located at (x, ?/) = (1, 3) 
and (2,3). 

The decomposition of py into two subspaces occurs due 
to rapid sign change of the magnetic moment at around 
the small island. The density matrix py is constructed 
by inner product between two column vectors of m, and 
then the inner product between the vectors for small and 
large islands vanishes due the sign change. Thus, the 
inner product means spatial correlation between the is- 
lands. As for the 3x3 submatrix, the absolute values 
of the off-diagonal components are as large as the diago- 
nal components. This is also clear, since the off-diagonal 
component comes from inner product between two col- 
umn vectors inside of the large island. Then, the eigen- 
values split into large one and zero, leading to Ai = 6. 
This situation is similar to that of the energy-level split- 
ting into bonding and anti-bonding states after mixing of 
two orbit als. 

By solving the eigenvalue equations, we obtain 



1 

71 



(39) 



and 



/ 1 



V 



,V2 



(40) 



Finally, each layer of m = m'^^^ is reconstructed by 
„^W = V6U,®%=(^1 ; :} 0), (41) 



m 



= V2U2^% 



1 
-1 



(42) 



Therefore, we find that the small and large islands are 
automatically decoupled into different layers. Realistic 
cases are more comphcated, but the feature of the rapid 
sign change still remains even in those cases. 



V. THREE-STATE POTTS MODEL 

We would like to know whether our results are gen- 
eralized for much broader universality classes. Here, we 
examine the 3-state {q = 3) Potts model: 



H 



(43) 



with dj = -1, 0, 1, (5(cr, a') = l{a = a'), -l{a 7^ a'). The 
central charge for Zq symmetric CFT and the critical 
temperature are c = 2{q — l)/{q-\-2) = 4/5 and Tc/J = 
2/log(l + y^) = 1.98994. We perform MC simulation 
to obtain snapshots, where we take maximally 10^ MC 
steps. We remember that we need to take at least L > 64 
to catch the critical behavior in the Ising model. 
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FIG. 12: (a)(b) T and L-dependence of 5' and a = 5' - logL 
for 3-state Potts model: L = 81 = 3^ (blue), L = 128 = H 
(purple), and L = 256 = 2^ (red). The data for L = 81 and 
L = 128 are avaraged by 10^ samples. We do not take the 
avarage for L — 256. 

Figures [121 (a) and (b) show T and L-dependence of 
S and OL = S — logL. The abrupt decrease in S occurs 
at around T = 1.99 J, suggesting that Tc is very close to 
the exact value. In Fig. [12] (b), we again find the scaling 
Eq. ([15]), 5' = logL + a, for T > T^, and ot = -7r/4 in the 
large-T limit and a ~ —2 at Tc. Therefore, Eq. ([T5]) is 
hold in the 3-state Potts model. The slope of the finite-L 
scaling for the a value at Tc is also different from that 
above Tc. The slope at Tc decreases with T, while the 
slope above Tc increase as L. The scaling at Tc, 5* = 
log!/ — 2, would come from critical behavior as we have 
already discussed in previous sections. 




FIG. 13: Temperature dependence of entanglement spectrum 
En. T = 2.0J (red), T = 1.99J - Tc (purple), and T = 1.98J 
(blue). We plot the data of 10 samples on the samle pannel. 
(a) L = 81 = 3^. (b) L = 128 = 2^ 

Figure [13] shows T-dependence of T^n for T = 81 = 3^ 
and T = 128 = 2^ near Tc. We find that the kink starts 
to appear with decreasing T, and finally the gap opens 



below Tc. Near the kink positions, ferromagnetic islands 
with two different length scales compete with each other, 
and then the data are somehow scattered. Thus, we use 
this scattering as a sign of the kink. The kink positions 
are located at nj^inj^ = 3 and Tikink 

= 8 - 9 for T = 

81 and L = 128, respectively. The result suggests that 
^kink loggT, and actually the logT dependence on a 
at Tc is related to critical behavior. 

The data at T = 2.0 J slightly above Tc also bend at 
n larger than Ukink- This is because T is one of energy 
scales, and thus is related to the length scale. Higher 
temperature represents shorter wave length, and thus the 
bending position shifts to larger n region as we increase 
T. Actually, in Fig. [13] (b), a weak kink is observed at 
^/cin/c = 14 15. We expect that this kink becomes clear 
as we increase L. 




FIG. 14: 5'x at T = 1.99 J - Tc for L = 256 (open purple 
circles) and L — 128 (filled purple circles). For comparison, 
5'x at T = 1.98 J is also plotted (blue circles). The slopes of 
the guide lines are 1/6 and 1. 

We finally show finite- x scaling for in Fig. [14] We 
could not conclude definitely the presence of the scaling 
relations (|27|) and (|28|) , but this size dependence looks 
similar to that in the Ising model. Therefore, we believe 
that the data asymptotically approach the scaling rela- 
tions previously obtained, and the scaling, S = log T — 2, 
comes from = log x — 2 for x ^ L in the large-T limit. 

VI. DISCUSSION AND SUMMARY 

Up to now, we have observed two scahng relations con- 
cerned with length scales of ferromagnetic islands in the 
snapshot at Tc. The first one is given by 

^x = ^logX + 7', (44) 

for X ^ ^/cin/c and a small positive 7' value. This scaling 
directly represents long-range part of critical behavior. 
The second one is given by 

= log X + 7 = ^ log ^kink + 7 + ' (45) 
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for X > "^kink- This leads to 

S = \ogL-2, (46) 

for X = ^ ^ CO. The factor /S.S^ represents correction 
to short-range spin correlation. Although Eqs. (|44|) and 
(|45]) are only numerically confirmed in the Ising model, 
Eq. (|46|) is hold in the Ising and the 3-state Potts model. 
Thus, their equations are expected to be fundamental 
ones irrespective of the value of the central charge. 

Equation (|44|) is quite analogous to the MPS scaling in 
Eq. (HI, but the definitions of x ^i"^ different with each 
other. In our case, x is the truncation number of the SVD 
for the snapshot of 2D classical spin configuration, while 
X represents the matrix dimension in the MPS. Although 
the MPS is also originated from the SVD, the role of x 
on the scaling is quite different. In the MPS scaling, we 
need more and more matrix dimensions for large-^ cases. 
On the other hand, in the present case, it is enough to 
take X = ^kink to catch the long-range correlation of the 
model. The relation between the length scale and x ^tre 
reversed in the MPS and our analysis. 

Usually, ID quantum sysytems can be mapped onto 
2D classical ones by the Suzuki- Trotter decomposition. 
We consider that our observation potentially catches ID 
quantum physics with the same universality class. Let 
us remember the transverse-field Ising chain that can be 
mapped onto the anisotropic 2D classical Ising model. 
Although we do not know the unique quantum corre- 
spondence of the isotropic Ising model, but we expect 
the presence of the quantum model. Then, multiplying 
a factor (l/Q)cA into S in Eq. (|46|) at Tc leads to 

gclassical ^ ^^J^g, (47) 

6 

with = L and this may correspond to the entanglement 
entropy of ID quantum systems in Eq. ([3]) (or Eq. (|2j)). 

Finally, we briefiy touch on the presence of the topo- 
logical term in our formalism [13, [HI . Equation (|23|) can 
be transformed into 

(^-^ log 2^ Ukink = ^ log L - log V2. (48) 

This equation may be composed of the deformed CFT 
scaling and the topological term. Namely, 

rikink (x^cAlogL - log ^/q, (49) 

where we assume A = 1. We start from the torus ge- 
ometry due to the periodic boundary condition, and cut 
the trus in order to calculate the reduced density matrix. 
Then, the number of the cut for each direction is 1. The 
relation (|^9|) seems to be consistent with the numerical 
data for the 3-state Potts model. When we take c = 4/5 
and (7 = 3, Eq. (|^9]) is given by Ukink ^ 41og3L — 15. 
Then, we obtain 41og3L — 15 = 1, 2.64, and 5 for 



L = 81 = 3^, 128 = 2^, and 243 = 3^ respectively. 
On the other hand, the kink positions in Figs. [13] (a) and 
(b) are located at 3, 8 ~ 9, and 14 ~ 15. Thus, we find 

^rikink ^ 4 log3 L - 15. (50) 

Therefore, Eq. (|49]) is satisfied for two models with the 
different magnitudes of the central charge, respectively. 

The second term in Eq. (|49|) . — logy^, would be 
a sign of topological nature in a ID quantum system 
corresponding to the present 2D classical Ising model, 
although usually the topological entanglement entropy 
Stopo is defined on the 2D quantum systems as 5 = 
aL + Stopo- However, possible presence of topological 
effects even in ID has been discussed recently in the Hal- 
dane phase of 6* = 1 chains [19]. According to the def- 
inition of the entanglement gap, the gap separates low- 
energy topological levels from high-energy generic ones. 
Therefore, the doubly-degenerate edge modes are topo- 
logically protected, and this protection is coming from a 
set of symmetries. In the present case, on the other hand, 
the lowest-energy state below the entanglement gap rep- 
resents ferromagnetic order. Then, the Zq symmetry is 
spontaneously broken. It would be important to clarify 
whether such a dual nature really exists or not, since the 
Haldane phase is not critical. 

In summary, we have examined the entanglement en- 
tropy of the coarse-grained MC snapshots of 2D square- 
lattice Ising and 3-state Potts models. Up to now, this 
type of entropy has not been examined yet, but we have 
found rich physical aspects hidden in the entropy. In 
particular near Tc, the entropy naturally produces the 
CFT scaling and the scaling for the entanglement sup- 
port of the MPS approximation on a unified framework. 
In addition to the entropy, the entanglement spectrum 
also gives us the scaling relations due to the presence of 
the multiple entanglement gap. A key ingredient of these 
scaling relations is fractal nature of ferromagnetic islands 
near Tc which have various length scale. Then, the SVD 
automatically decomposes the original snapshot into a 
set of images with different length scales, respectively. 
The kink structure of the entanglement gap character- 
izes this decomposition, and further suggests the possible 
presence of the topological term. Based on the present 
results, we are very interested in the detailed understand- 
ing of the duality between ID quantum and 2D classical 
systems. Examinations of the scaling and the duality in 
much broader universality classes will be important fu- 
ture works. 
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